function[cons]=invutility(beta,rho,U);
if rho==1
cons=exp((1-beta)*U);
else
    cons=((1-beta)*(1-rho)*U+1).^(1/(1-rho));
end